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ABSTRACT 
A  theory  of  low-level  turbulence  in  compressible  plasmas  is 
presented.  It  is  characterized  by  a  separation  of  scales  between  the 
mean  fields  and  the  fluctuations,  which  are  of  low  amplitude.  Use  of 
conservation  laws  of  the  system  enables  the  derivation  of  important 
factors  in  the  functional  form  of  the  turbulent  diffusivity 
coefficients.  In  particular,  it  is  shown  that  turbulent  heat 
conductivity,  like  its  molecular  counterpart,  transports  heat  in  the 
direction  opposite  to  the  mean  temperature  gradient.  The  turbulent 
resistivity,  on  the  other  hand,  acts  non-classically  and  causes 
flattening  of  the  parallel  current  profile  so  that  the  plasma  tends  to 
approach  a  "Taylor  state".  The  case  of  planar,  viscous,  fluid  flow  is 
treated  in  detail  as  a  special  example.  It  is  shown  that  turbulent 
viscosity  in  the  low-level  turbulence  case  simply  causes  enhancement  of 
the  molecular  viscosity. 


•Presented  at  the  Fifth  International  Beer-Sheva  Seminar  on  MHD 
Flows  and  Turbulence,  March  2-6,  1987,  Jerusalem,  Israel.  To  appear  in 
the  AIAA  Progress  Series. 
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I.   INTRODUCTION 

The  subject  of  turbulence,  namely,  the  description  of  flows  which 
vary  rapidly  in  space  and  time,  does  not  easily  lend  itself  to  either 
analytical  or  numerical  methods.  Indeed,  many  of  the  statistical 
methods  •  used  by  physicists  and  engineers  are  not  established  on  a 
rigorous  basis  and  are,  in  fact,  known  to  yield  incorrect  answers  under 
some  circumstances  which  are,  unfortunately,  not  always  definable. 
There  is  an  exception  to  this  rule  and,  quite  remarkably,  it  has 
escaped  attention  until  very  recently.  This  is  the  area  of  low-level 
turbulence. 

By  the  term  "low-level  turbulence"  we  mean  that  the  ratio  of 
amplitudes  of  the  fluctuations  to  the  mean  flow  is  small,  as  well  as 
that  there  is  a  clear  separation  of  both  temporal  and  spatial  scales 
between  these  two  components  of  the  flow:  the  fluctuations  vary  on  a 
faster  time  scale  and  a  shorter  spatial  scale  than  the  mean  flow. 
These  assumptions  make  possible  the  use  of  multiscale  asymptotic 
methods  which  lead  to  simplified  equations  and  allow,  in  some  cases, 
the  derivation  of  an  exact  expression  for  the  turbulent  "eddy" 
diffusivity.  An  example  is  given  in  Section  II,  In  this  context  we 
mention  the  work  of  McLaughlin  et  al.-^  who  use  such  methods  mostly  for 
fluctuations  of  large  amplitude.  This  limits  the  extent  of  the 
analytical  results  that  can  be  derived. 

Low-level  turbulence  may  not  be  too  pressing  a  subject  in  fluid 
dynamics,  but  it  is  precisely  the  natural  turbulence  problem  in  plasma 
confinement  studies.  Experience  shows  that  the  confinement  of  a  hot 
plasma  by  a  magnetic  field  is  a  very  delicate  matter.  If  the  plasma 
undergoes  large  fluctuations,  confinement  is  lost  roughly  within  a 
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fluctuatlon  period.   By  definition,  therefore,  fluctuations  in  a 
confined  plasma  are  of  small  amplitude.  The  separation  of  scales  too 
arises  naturally  when  one  identifies  the  origin  of  the  fluctuations  in 
the  myriad  of  instabilities  affecting  the  plasma,  many  of  which  tend  to 

be  localized  about  a  "resonant"  magnetic  field  line.  "^   Perhaps  this 

6— R 
recognition  gave  rise  to  some  works     which  take  advantage  of 

disparate  scales   in  the  plasma  dynamics.  Their  results  are  still  too 

cumbersome  to  offer  an  understanding  of  the  turbulent  dlffusivity  of 

the  plasma. 

A  most  intriguing  phenomenon  which,  we  believe,  is  explained  by 

the  effect  of  low-level  turbulence  is  the  sfxjntaneous  reversal  of  the 

toroidal  magnetic  field  in  Reversed-Field  Pinch  (RFP)  experiments.' 

According  to  the  theory  of  plasma  relaxation  due  to  Taylor    (and 

recently  reviewed  by  him   ),  a  slightly  resistive  plasma  tends  to  reach 

an  equilibrium  state  of  minimum  magnetic  energy  subject  to  the 

conservation  of  its  magnetic  hellcity,   thus  producing  a  force-free 

state,  now  known  as  a  "Taylor  state",  with  VxB  =  X^B,  where  B  is  the 

magnetic  field  and  X^  is  a  constant.   Solutions  of  this  equation 

approximate  general  features  of  observed  magnetic  configurations  in  the 

RFP,  including  toroidal  field  reversal  ,  but  the  theory  does  not 

provide  a  dynamical  description  of  the  relaxation  process.  A  numerical 

calculation  due  to  Strauss,   and  based  on  a  magnetohydrodynamic  (MHD) 

description  of  the  plasma  with  disparate  scales  taken  into  account, 

does  show  such  an  effect.  Our  analytical  work  with  Bhattacharjee  on 

low-level  turbulence  in  a  plasma,  -"         successfully  determined  an 

explicit  expression  for  the  eddy  resistivity  (which  involves  a  fourth 

derivative  of  B),  along  with  its  dependence  on  the  mean  fields. 


Varlous  other  considerations^"^ ^  lead  to  similar,  If  less  precise, 
results.  It  was  demonstrated^'   that  a  plasma  under  the  Influence  of 

such  eddy  resistivity  does  tend  to  approach  a  Taylor  state.  It  was 

1 B 
also  shown   that  the  plasma  need  not  obey  Taylor's  energy  principle. 

If  anything,  it  follows  some  entropy  production  principle. 

In   this   paper  we  extend  our  earlier  work,^^'    in  which 

resistivity  was  considered  to  be  the  only  diffusive  effect,  to  plasma 

with  finite  resistivity  and  heat  conductivity.  This  also  means  that, 

as  in  Ref.  1^4,  we  deal  here  with  a  compressible  plasma,  a  case  which 

is  almost  never  treated  in  the  literature  of  the  analogous  question  of 

1 Q 
the  "dynamo"  problem.  ^  Rather  than  getting  an  exact  expression  for  the 

"eddy"  heat  conductivity,  we  are  content  here  with  producing  an 

expression  for  it  in  terms  of  means  of  fluctuating  quantities,  and  with 

exploring  its  general  structure.   It  turns  out  to  be  more  similar  to 

the  molecular  heat  conductivity  than  eddy  resistivity  Is  to  its 

molecular  counterpart.   This  Is  done  in  Sections  IV  and  V,  after  a 

review  of  earlier  results  in  Section  III.  But  first  we  start  with  a 

demonstration  the  multlscale  method  and  its  effectiveness  in  an  easier 

problem:  two-dimensional  turbulent  flow  and  the  question  of  eddy 

viscosity.   In  this  special  case,  an  expression  for  the  eddy  viscosity 

can  be  explicitly  determined  if  the  level  of  turbulence  is  low. 


-5- 

II.   TWO-DIMENSIONAL  FLUID  FLOW 

We  consider  the  planar  motion  of  an  Incompressible  fluid,   which 
obeys  the  Navier-Stokes  equations 

3v         , 

_1  ♦  vVv  +  _  Vp  -  vAv  (1) 

3t   -  -   p 

V-y  -  0  .  (2) 

Here  v,  p  and  p  are  the  velocity,  pressure  and  the  constant  density, 
respectively,  and  v  is  the  constant  viscosity.  The  curl  of  Eq.  (1) 
yields 


—  +  vVii)  ■=  vAto  ,  (3) 

3t 


where  u  is  the  component  of  the  vorticity  Vxy  normal  to  the  plane. 

Let  us  write  every  quantity  as  the  sum  of  a  mean  part,  denoted  by 
a  subscript  o,  and  a  fluctuating  part,  denoted  by  a  subscript  1.  Thus, 
V  =  v  ♦  v^  ,  where  v^  =  <v>  and  <v^>  ■=  0,  and  the  brackets  <  >  denote 
either  an  ensemble  average  or  an  average  over  the  spatial  and  temporal 
scales  of  the  fluctuations.  We  now  make  the  following  assumptions:  (1) 
the  mean  state  is  near  steady  state  and  evolves  slowly  on  a  viscous 
time  scale,  while  it  varies  spatially  on  the  scale  a  of  the  system, 
i.e.,  aly^l  -  0(1),  a^3/at  =  0(v),  a\V\  -  0(1);  (2)  the  fluctuations 
vary  on  a  shorter  space-scale  and  a  faster  time-scale,  i.e., 
a^3/3t  >>  0(v),  a|V|  >>  0(1);  (3)  the  fluctuation  level  is  so  low  that 
|(i)J  «  IujqI  and  Vu^  is  at  most  order  1,  i.e.,  a|VwJ/wQ  -  0(1).  In 
the  following  we  will  consider  all   quantities   to   have   been 
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nondimensionalized;  spatial  variables  scaled  with  a,  and  velocities 
scaled  with  U,  where  U  is  a  typical  mean  velocity. 
Averaging  Eq.   (3)  yields 


3(1) 

-r—   *  yo'^"o  *  f  "  vAiOq  ,      F  =  div  <Li■^■V'^>    .  W 


The  term  F  embodies  the  effect  of  the  fluctuations  on  the  mean  flow 
and,  according  to  our  ordering,  should  be  small.  To  leading  order  the 
mean  flow  is  in  steady  state,  satisfying  Vo'^^^  -  0.  Thus,  the 
evolution  may  be  viewed  as  taking  place  through  a  sequence  of 

near-equilibrium  states.     This  is  analogous  to  the  Grad-Hogan  theory  of 

20 
slow  plasma  diffusion.    For  simplicity  let  us  assume  that  each  such 

steady  state  generates  a  nested  family  of  streamlines  as  in  Fig.   1. 

To  study  F,  let  us  write  Eq.  i^)    in  the  form 

A  (^  ^l)    .   div  (oi^v^  -  vu)^V.3^)  -  ..^F  =  -v|Va,j2  .     (5) 

We  can  generate  another  such  equation  by  writing  Eq.   (3)  as 

—  (1  u)^]  +  div  (I  uj^v  -  vwVu)]  =  -vIVwl^  ,  (6) 

it     2  2       ~  '   ' 

and  averaging.  The  result,  to  order  v,  is 

A   (I  0,^)    *   div    (1  oj^Vq  -   vajoVu)^]    *   div   <a)oa)iyi    *  1  (ofyQ> 

-  -vIVuqI^  -  v<|Vu)J^>    .  (7) 
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Note  that  the  averaging  Is  taken  under  the  derivatives,  and  that  the 
term  3<(i)^>/3t  is  o(v).  Comparing  Eqs.   (5)  and  (7),  we  have  to  0(v) 


uJqF  -  v<|Va)J^>  ♦  div  <u)Qa)^y^  ♦  _  u^\^>    .  (8) 


Writing  F  -  V'5,  where  C  "  *^'*'^y^>,  we  have  (if  there  is  no  flow  across 

the  boundary)  /Fd  x  -  0,  while  from  (8)  we  have  /tu^Fd  x  >^  0.  The 

integrals  here  are  taken  over  the  flow  region.  Integrating  by  parts, 
we  get 


j  (§-Vu)Q)d^x  <  0  .  (9) 


20  ?1 
In  analogy  with  the  Grad-Hogan  theory,   '   we  consider  the  geometry  of 

closed  streamlines  of  the  mean  flow,  label  each  curve  by  A,  the  area  it 

encloses,   A  =  A(x),   and  write  d  x  -  dAd!L/|VA|,  where  dX,  is  a  length 

element  along  the  closed  curve.  Thus,  Eq.   (9)  can  be  rewritten  as 


d(i). 


\''-wK "--''-'  •  ^''^ 


Here  E,^  is  the  component  normal  to  the  curve,  the  subscript  A  indicates 
the  particular  curve  with  label  A,  and  we  have  used  the  fact  that,  to 
0(1),  L^^  «=  oi^CA).  Now,  if  the  correlation  of  the  turbulence  across 
mean  streamlines  is  negligible,  as  is  the  case  when  the  multiscale 
method  applies,  then  inequality  (10)  must  be  satisfied  separately  on 
each  curve.  We  thus  find  the  following  property  of  the  eddy  viscosity 
term: 
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tn^^--K^T^.  (11) 


Ja  "        dA  ' 

where  K  is  a  non-negative  function.  This  conclusion  will  shortly  be 
verified  in  the  example  given  below.  It  is  worthwhile  to  point  out 
that  similar  arguments  could  have  led  to  the  conclusion  that  iFd£/|VA| 
should  be  proportional  to  Wq,  or  to  other  conclusions  along  similar 
lines.  We  emphasize  that  conclusion  (11)  relates  to  §  itself,  and  not 
to  its  derivatives.  For  turbulence  with  short  correlation  length, 
derivatives  of  §  can  be  much  larger  than  §  and  an  ordering  scheme 
leading  to  Eq.   (8)  might  fail. 

The  implication  of  (11)  can  be  seen  if  we  take  the 
(microcanonical)  average^°'^^  of  Eq.  (J4),  that  is,  integrate  it 
d2./|VA|.   Using  incompressibility,  and  the  fact  that  for  any  f 


where  the  subscript  A  on  the  right-hand-side  indicates  integration 
within  the  curve  A,  we  find 

^"c'«.t)«^|,  t„»t.v^|^|   |,.|.t|.  (.2, 


Where  the  first  term  follows  from  non-trivial  averaging.  Equation 
(12)  is  a  one-dimensional  diffusion  equation  for  ^^(A.t),  with  respect 
to  the  given  geometry  of  streamlines.    (Of  course,  the  streamlines 
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themselves  also  change  at  the  same  time,  and  the  slow  evolution  problem 
needs  to  be  solved  Iteratively.   •   )  Combining  Eqs.   (11)  and  (12),  we 
get 

±.o(A.t)-A|[v|jVA|d.  .K^]^}  .  (13) 

2 
The  K   term,  which  may  depend  on  the  shape  of  the  streamlines, 

expresses  the  enhancement  of  diffusion  by  eddy  viscosity.  Its  effect, 

In  this  case,  appears  to  be  similar  to  that  of  the  molecular  viscosity. 

Using  similar  arguments,  we  will  see  in  the  next  section  that  the  same 

Is  not  true  for  the  electrical  resistivity. 

It  will  be  comforting  to  derive  property  (11)  directly  from  the 

dynamical  equations  of  the  fluctuations.  For  this  we  need  to  adopt  a 

more  specific  ordering   (which  is  why  we  preferred  to  derive  the  more 

general  result  first).  We  use  the  simplest  geometry  possible  and 

assume  the  mean  streamlines  to  form,   to  leading  order,  concentric 

circles.   In  polar  coordinates  r,e,  we   have   Uq  ■=  cOq^'"^  *  0(v). 

Subtracting  Eq.    C*)  from   (3),  we  get  the  evolution  equation  of  the 

fluctuations. 


3u. 

— —  +  ^o'^'^i  *  -1  *'"o  *  -1  *^"i  ~  ^^^  '^'^'lyi^  "^  '^^'^1  •  ^^^^ 


We  assume  the  fluctuations  to  vary  rapidly  across  mean  streamlines,  but 
to  have  0(1)  variation  along  them.   Introducing  the  smallness  parameter 
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6,  we  use  the  small  scale  variables  p  -  r/6,  i   -  6  t,  and  make  the 
multiscale  ansatz 


V  -  y,(p,T;r,e,t)  ,  (15) 


where  we  take  v  -  0(6  ).  Note  that  in  Eq.  (I^O  we  now  have  to  replace 
3/9t  ■»  3/9t  +  6^9/91,  9/9r  ■►  9/9r  +  6~^9/9p.  The  order  of  the 
perturbed  velocity  components  (indicated  by  subscripts)  is  taken  to  be 


^1r  "  S^v^^)  +  0(6^  ,      v^g  -  6^v^^)  *    0(63) 
..,  -  6^(1)  .  6^2)  .   6^3)  .  0(6^  . 


where  w    -  9vg  /9p,   but  higher  order  terms  of  tu^  depend  on  higher 
order  terms  of  the  velocity.   Up  to  0(63),  £q^   (^i^)  reads 


.-3(^  •  .V-V.,)  .  ^41   .  .;v<3>  .  (v<3,  ^  .  -  .<.)  ^)„(-.  .  „  3^  .   „e, 


Here  w^  -  doj^/dr,  and  we  notice  that  the  term  div  <uiy^>  in  Eq.  (1^), 
which  is  of  the  order  of  w-jV.  •=  0(6  ),  has  dropped  out.  The  terms  in 
the  first  parentheses  in  Eq.  (16)  are  sums  of  terms  ranging  from  0(6) 
to  0(63).  Thus,  to  0(6),  we  have 


(^  +  w-V]J^^   -  0  ,  (17) 


•9t 


-o 
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and  the  same  holds  for  u     K       Equation  (17)  implies  that  w^  '  Is 
advected  circularly  by  the  mean  flow,  as  should  be  expected. 

It  is  useful  to  replace  6  by  the  angular  Lagranglan  coordinate  of 
the  mean  flow,  e^  -  6  -  fiQ(r)t,  where  AqCt)  is  the  rotational  frequency 
of  the  mean  flow  on  the  curve  r  -  const.  The  corresponding  change  in 
Eq.  (16)  is  3/36  -»  3/39^  and  3/3t  ♦  v^'V  -»  3/3t.  At  0(6^),  the 
structure  of  Eq.   (15)  is 


3J3) 


"sT"  '  ^^^•■^•'"•^o^  •  ^^S) 


where  the  independence  of  the  right-hand-side  on  t  follows  from  Eq. 
(17)  and  the  fact  that  Vg  ^  and  v^^^  are  determined  from  dj  '  by  using 
the  definition  of  u    and  the  divergence-free  condition. 


_il_*I_i_=0.      .J^^J^^.  (19) 

3p    r  36^  3p 


The  solution  w^^^  of  Eq.   (18)  grows  linearly  in  t,  unless  f  =  0.   We 
thus  require 


9"^^^  .   •..(^)  .  f..(^)  3  .  1  ..(2)   3  >  (1)     3^u)^^^ 


+ 


'v(3)  .  (v(3)  3  .  I  ,(2)   a  )J1)  „  ,  r^   _    (2Q^ 


Such  a  non-secularity  condition  is  common  in  multiscale  analyses. 

Solutions  of  Eqs.   (19), (20)  depend  on  the  independent  variables  t,   p 

» 

and  Oq,  as  well  as  on  the  three  parameters  v,  r  and  w^.  The  parameters 
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can  be  eliminated  by  the  following  scaling  transformation  (which 
effectively  sets  them  to  1): 


(21) 


•^o  ^"^o 


v^3)  .  [1.   v3.;)l/\(3)  ,     v(2)  .  (rvu>;)1/2,(2)  , 


The  eddy  viscosity  term  can  now  be  determined  to  order  v  -  0(6  ), 


<a),v,^>  -  -CvUq  ,     C  >  0  ,  (22) 

where  C  is  some  universal  constant.  Its  positivity  follows  from  Eq. 
(9)  or,  more  directly,  by  multiplying  Eq.  (20)  by  o)^  and  averaging 
over  the  Independent  variables.  In  this  simple  geometry,  Eq.  (M) 
implies  that  the  eddy  viscosity  simply  equals  Cv.  A  similar  relation, 
up  to  a  dependence  of  C  on  the  geometry  of  mean  streamlines,  is 
expected  in  the  more  general  two-dimensional  flows. 
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III.   TURBULENT  DIFFUSION  IN  INCOMPRESSIBLE  PLASMAS 

In  this  section  we  review  some  results^'  concerning  the  eddy 
resistivity  of  Incompressible  plasmas.  The  next  section  develops 
additional  results  for  the  compressible  case,  Including  an  expression 
for  the  eddy  heat  conductivity.  The  plasma  is  assumed  to  obey  the 
resistive  MHD  equations 

p[°    +  vV)v  +  Vp  -  JxB  (23) 

3t    ~    ~        —  — 

3B 

-1  +  Vx(nJ-vxB)  •=  0  (2^4) 

dt  -  -  - 

V«v  ■=  0  ,    p  ■=  const.  (25) 

where  n  is  a  "small"  constant  scalar  resistivity,  and  J  ^  VxB  Is  the 
current  density.  The  plasma  is  assumed  to  be  confined  and  isolated 
inside  a  toroidal  solid  conductor,  on  which  the  boundary  conditions 

B«n  -  0  ,    v-n  =  0  ,  (26) 

hold,  where  n  is  normal  to  the  boundary.   Similar  to  the  previous 

section,  we  consider  a  slowly  varying  mean  state,  but  this  time 

containing  only  diffusion-induced  flows,   upon  which  is  superimposed 

low-level  turbulence.  More  precisely,  and  using  previous  notations,  we 

use  the  following  assumptions:  (1)  for  the  mean  state,  ^\y.o\  "  ^^'^)» 
a^3/9t  =  0(n).   ajVl  =  0(1);   (2)  for  the  fluctuations,  a^3/3t  >>  0(n), 

alV|  >>  0(1);  and  (3)  the  energy  of  the  fluctuations,  be  it  magnetic, 

kinetic  or  thermal,   is  smaller  than  the  mean  magnetic  energy.  I.e., 


|bJ^,  p|yj^,  p^  <<  |Bq|^.  Since  J,  -  '«li.  the  fluctuating  current 
may  be  large.  We  allow  it  to  be  of  the  same  order  as  Jq.  In  the 
following,  we  will  consider  our  variables  to  have  been 
non-dimensional ized  with  respect  to  typical  mean  magnetic  field 
parameters. 

Averaging  Eqs.   (23)  and  (2^),  we  get 


^Po  ■  ^o-^o  •  ^27) 
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^  *  ^"(^Jo  -  Yo^lo  -  <}L}'h>^  '   0 


(28) 


where  we  note  that  Eq.   (27)  is  only  correct  to  0(1),  and  implies  that 
Bq  is  in  static  equilibrium.  Equation  (28)  then  requires  the  evolution 

to  take  place  through  a  sequence  of  equilibrium  states.   According  to 

20 
the  Grad-Hogan  theory,    the  velocity  v^  is  determined  (actually, 

adjusts  itself)  so  that  this  requirement  is  satisfied.   It  will  be 

assumed  in  the  following  that  B^  forms  a  family  of  nested  toroidal 

magnetic  flux  surfaces. 

We  now  derive  some  properties  of  the  turbulence-induced  electric 

field,  e  ^  ^^i^Bi^-  The  following  relations  must  hold  to  order  ri. 


£-Jo^^  - 


n<jf>dT  ,     I  e-B^dT  -  0  ,  (29) 


where  the  integrals  are  taken  over  the  plasma  volume.  The  first 
relation  may  be  shown  by  averaging  the  energy  equation 


-15- 


•^  (-^  B^  *  -^  Py^}  *  V-[^  Py^y  *  py  ■►  (TiJ-yxB)xB]  +  TvJ^  -  0  ,   (30) 

which  may  be  obtained  from  Eqs.   (23)-(25).  The  average  of  Eq.   (30) 
yields,  to  0(n), 


■^  ^^  §0^  *  ^o  *  '^H^^>   *  ^'£1  -  0  .  (31) 


where  c-j  represents  energy  flux.   Another  equation  for  the  mean 
magnetic  energy  may  be  obtained  from  Eq.   (28).  The  result  is 


^  (^  ?§)  *  ^Jo  -  S'io   *  ^-^2  ■  0  (32) 


The  surface  integral  /(c^-C2)*ndS  represents  the  energy  flux  due  to  the 
fluctuations  which,  for  an  isolated  plasma,  vanishes.  A  comparison  of 
Eqs.   (31)  and  (32)  yields  the  desired  property  of  e^J^. 

The  second  property  in  Eq.    (29)  may  be  proved  similarly  by 
averaging  the  helicity  equation  (with  B  -  VxA) 


3  r  9A   ^ 

-I.  (A'E)  +  2nJ'B  +  V«[2(tiJ-vxB)«A  +  _I  kA]  -  0  .         (33) 
3t   -  -      -  -     *■----   8t   - 


Averaging  yields,  to  0(n), 


■^   (A^-?o)  *  2nJo-?o  *   '-£3  =  0  ,  (3^^) 
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where  c^  represents  heliclty  flux.  However,   from  Eq.    (28)  we  can 
write 


dt 


■^   (Ao*?o^  *  2nJo-Bo  -  Se-Bq  +  V-Ci^  -  0  .        (35) 


J{cj-c^)-n6S     represents  the  heliclty  flux  due  to  the  fluctuations, 
which  again  vanishes.  Comparing  Eqs.   (3^)  and  (35)  yields  the  form 


S'§o  '  ^'S  •  ^^^^ 


where  t   has  zero  net  flux  across  the  boundary,  and  the  desired  result 
follows. 

From  the  two  properties  expressed  in  Eq.  (29)  we  can  guess  the 
form  of  the  component  of  e  parallel  to  B^,  along  the  same  lines  used  in 
the  previous  section.  For  a  low-pressure  plasma  (which  is  the  usual 
experimental  situation),  Eq.  (27)  implies  that  J^  is  nearly  parallel 
to  Bq.  Writing  Jq  -  XB^  and  using  Eq.  (36)  and  the  first  relation  of 
Eq.   (29),  we  get  after  integration  by  parts 


h 


C'VXdT  >  0  .  (37) 


This  is  analogous  to  Eq.    (9)  and  similarly  suggests,   if   the 
correlation  of  fluctuations  across  flux  surfaces  is  negligible,  that 


C,dS  .  k2  £i  ,  (38) 
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where  here  the  volume  V  labels  flux  surfaces  and,  again,  K^  >^  0. 
Notice  that  V«Jq  -  0  Implies  that  X   Is  nearly  constant  on  flux 
surfaces. 

Relation  (38)  was  derived  In  Refs.  13  and  ^U  directly  from  the 
dynamical  equations  of  resistive  tearing  modes.  It  is  Important  to 
notice  that,  when  substituted  in  Eq.  (28),  it  enters  as  a  diffusive 
term  of  fourth  order  in  spatial  derivatives  of  B^  (since  \,  the 
parallel  current,  is  itself  related  to  derivatives  of  B^).  In  this 
respect,  it  is  significantly  different  from  the  eddy  viscosity  term 
derived  before.   Indeed,   it  can  explain  the  phenomenon  of  plasma 

relaxation  to  a  near  Taylor  state,   '   which  cannot  be  explained  by  a 

22 

simple  enhancement  of  classical  resistivity. 

As  an  illustration,  we  present  here  a  numerical  calculation  of  a 
steady  state  solution  of  a  cylindrical  plasma  under  the  influence  of 
low  level  turbulence  and  an  externally  applied  electric  field  in  the  z 
direction  (which  is  the  plasma  axis).  We  solve 

nJ  -  E  -  Ez  ,  (39) 


where  the  constant  E  is  the  applied  axial  voltage  per  unit  length,  and 
c  has  the  form 


e  -  a^  -J-  (na^rB^x'  )'b  +  e^  -^  r^B  .  (^0) 

rB"^  B"^ 


Here  all  quantities  are  mean  field  quantities,  primes  denote  d/dr,  a  is 
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the  plasma  radius,  a^  and  Bq  are  constants,  and  the  pressure  p  is 
obtained  from  the  equilibrium  relation 


Or      2       r 


with  p(a)  -  0.   The  boundary  condition  dX/dr  -  0  at  r  -  a  is  used  to 
conform  with  properties  (29). 

Expression  (.HO)   was  derived  in  Ref.   1^  under  the  assumption  that 

h 
tearing  modes   are  mainly  responsible  for  the  component  of  e  parallel 

ii 
to  B,   and  interchange  modes   are  mainly   responsible   for   the 

perpendicular  component.  The  expression  given  here  for  the  latter  is, 

in  fact,  new  and  corresponds  to  a  version  of  the  mode  dynamics 

appropriate  for  low-pressure  plasmas.  -^  Interestingly,  this  component 

of  £  equals  -BqH  times  the  perpendicular  component  of  J.   Thus,   it 

corresponds  to  a  simple  enhancement  of  the  perpendicular  resistivity. 

Equations  (39)-(.H^)     are  solved  after  non-dimensionalization.  We 

define  the  volume-averaged  B  to  be 


B,  -  4:  f^  rB^dr  ,  (^2) 


and  scale  B  by  B^,  p  by  B^,   and  E  by  B^n/a.   In  addition,   r  is 

normalized  to  a.   Effectively,  we  now  have  a-  1,  n-  1,  0_<r£1, 

with  boundary  conditions  X  (1)  -  0  and  B^  -  1 .   In  the  calculation 

presented  in  Fig.   2,  we  take  a^  -  1,  as  has  been  estimated  for  RFP 

1  U 

plasmas,   and  E  -  12,  which  also  corresponds  to  typical  machine 
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1  U 
parameters.    6^,  here  taken  to  equal  100,  was  chosen  large  enough  to 

yield  a  sufficiently  low  plasma  pressure.  As  to  the  results,  Fig.   2a 

shows  the  celebrated  field  reversal,  which  is  not  possible  If  a^  -  0, 

that  is,  without  effect  of  the  parallel  eddy  resistivity.   Figure  2b 

demonstrates  that  the  equilibrium  is  indeed  approximately  a  Taylor 

state,  with  a  very  flat  X(r)  profile. 
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IV.   TURBULENT  DIFFUSION  IN  COMPRESSIBLE  PLASMAS 

Extension  of  the  previous  analysis  to  the  Ceise  of  a  compressible 
plasma  is  not  difficult,  but  it  requires  a  redefinition  of  e.  Such  a 
device  has  already  been  used  in  Ref.  ^U.  However,  the  guiding 
consideration  given  there  is  inappropriate  if  the  plasma  has 
non-vanishing  diffusion  coefficients  other  than  the  resistivity,  a  case 
not  treated  in  Ref.  1^.  This  situation  is  corrected  in  the  present 
section,  and  we  are  able  to  deal  with  a  compressible  plasma  subject  to 
the  effects  of  heat  conduction  and  viscosity,  as  well  as  resistivity. 
The  plasma  evolves  according  to  the  following  equations: 


p(_  *   vV)v  ♦  Vp  -  JxB  *  V'ivo)  CiB) 

3t   ~   ~        -  ~       « 

8B 

—  +  Vx(nJ  -  vxB)  =  0  (^U) 
3t       -   -  - 

l£  *   V.(pv)  =  0  (^5) 

St 

—  (ps)  *  V.(psv)  -  III  u  tr  (o^)  *  111  nJ-J  *    (>-1)V-(l  kVt)  *  ^   kVT-VT  ,   (iJ6) 

3t  -      21         m  I  ~    ~  T  jt 


where  the  specific  entropy  s  satisfies  the  equation  of  state 
s  •=  log  (p/p^).  Here  T  =  p/p  is  the  temperature,  "T  is  a  constant,  o  is 
the  rate-of-strain  tensor  with  components 


3v.   3v.   p 
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and  tr  Indicates  the  trace  of  the  tensor.  Note  that  the  last  two  terms 
of  Eq.   (^46)  are  written  In  a  way  which  enables  easy  averaging. 

We  consider  here  the  case  of  a  tensor  resistivity  ti  (assumed  to  be 
small),  and  a  tensor  heat  conductivity  ic,  while,  for  simplicity,  the 
viscosity  u  is  taken  to  be  a  scalar  function.  All  the  diffusion 
coefficients  are  allowed  to  depend  on  B,  p,  and  T,  but  not  on  their 
derivatives  (which  is  the  physical  situation).  We  assume  now  that  J, 
and  VT^  are  at  most  0(1),  that  tc  is  of  the  same  order  as  n,  and  that  y 
is  so  small  that  vi|Vv.|j^|  is  at  most  0(ti),  for  any  component  i  of  the 
fluctuating  velocity.  These  assumptions  imply  that  the  fluctuating 
parts  of  n,  <  and  y  are  too  small  to  enter  the  averaged  equation. 
Thus,  whenever  the  diffusion  coefficients  are  used  in  the  following, 
they  should  be  taken  as  evaluated  for  the  mean  fields.  We  follow  again 
the  general  ordering  scheme  spelled  out  in  the  previous  section,  with 
the  added  requirement  that  |pi|  <<  p^,  except  that  y^  is  allowed  to  be 
as  large  as  v^ . 

While  averaging  equations,  it  is  desirable  to  average  them  in 
their  conservation  form  for  then  the  averaging  operation  is  carried  out 
under  the  derivatives  and  the  order  of  the  various  terms  is 
transparent.  The  average  of  the  momentum  conservation  equation,  Eq. 
(i)3)  combined  with  (A5),  yields  again  to  0(1) 

^Po  ■  ^o-lo  •  ^^7) 

while  the  averages  of  Eqs.   (iJH)-(i)6)  yield  to  0(n), 
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-_  +  Vx(nJo  -  yo«Bo  -  <v^xB,>)  -  0  ,  (iJ8) 

-^  *   V-Cp^Vq  ♦  <PlYi>)  -  0  ,  (i<9) 


±  (Po^o)  *   ^-(Po^oifo  *  Po<Siyi>  *  So<PiVi>) 


111  u  tr  <o^  *  111  (nJo'Jo  +  <riJT-J,>) 


*  (T-1)V.[_  kVIq)  ♦  12.   UVT^'VT^   *  <.cVT,-VT,>)  ,     (50) 
o        Tr 


where  <p^s^v>  was  dropped  from  Eq.   (50)   as  too  small.   Using  Eqs. 
('43)-(^6)  we  form  the  energy  equation 


9t  '2  -    2  -    Y-1 


+  V«[(tiJ-vxB)xB  *  -  py^v  ♦  py  -  kVT  -  voy]  «  0  .    (51) 

^  T""  I  K*" 


Averaging  it  yields  to  0(n) 


^  ^^  -0  *  TT^  *  ^-i^^^-o  -  Yo'-lo  -  <Yl''Bi>)xB^  *  1  Po<Y^> 


where  we  still  use  y  '  y^  +  y^ ,  and  a  term  of  order  Pi|y|   is  dropped 
as  too  small.   The  term  (y''B^)xB^  is  also  dropped  for  the  same  reason. 
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Comparing  Eq.   (52)  with  an  expression  for  3(Bq/2  ♦  PQ/(T-1)]/3t  based 
on  (K7)-(50),  we  find 


<v,«B,>-J^  *  T — L —  <PiVi>«7T^  -  —2-  <p,v->.Vs, 
_1  -1   -o   (Y-l  )T„  ^-'    o   Y-1    '-'    ' 


ii  tr  <of>  -  <nJ-,*Ji>  -  —  <tcVT^«VT,>  +  f  , 
2     .       "   ~     "^o 


(53) 


where  we  have   used   the   relation   (to   the   relevant   order) 
s^  -  P^/Po  -  ^Pi/Po'  S"'^  where 


f  =  -V.<(p^  *  !o*!l^^1  ■  (yr?1^?o  *  \   Po^V  •  (5iJ) 


The  function  f  appeared  in  Ref.  \^  ,  where  it  was  argued  that,  if  any, 
only  the  last  term  of  (5^)  which  expresses  kinetic  energy  flux,  may  be 
large  enough  to  enter  the  mean  equations.  The  quantity  p^  *  §o*§l  •  ^'^^ 
fluctuating  total  pressure,  is  smaller  than  B^  for  small-scale  motins, 
while  the  second  term  does  not  involve  flux  across  the  mean  magnetic 
flux  surfaces. 

We  can  simplify  the  mean  equations  considerably  by  using  the 
following  quantities  (which  will  be  discussed  at  the  end  of  this 
section) , 


E  -  <v,xBi>  -  J-  <PiVi>«B„  (55) 

-1  -1      p^    1-1    -o 


-2H- 
ilo  -  ^o  *  7-  <Piyi>  •  (56) 


The  evolution  of  mean  quantities  is  now  determined  by  Eq.    (^7),   as 
well  as  by 


3B 

-^   *  VxCnJ^  -  z  -   u^xB^)  -  0  (57) 

9Po 

-gT  *  ''^^0^0^  ■  0  (58) 

-g^  *  i^o'^Po  -^  ^Po^-i^o  ■  (^-l^t^Jo-io  -  i'io   *  '-(^VT^)  -  V-(^  <SiyT>)  *  f]   (59) 


The  relation  p^  ■=  p^T^  also  holds  to  0(1),  which  is  the  required  order 
for  Eqs.   (57)-(59).  Equation  (53)  now  takes  the  form 

i'io  *   ^  Po^^li^l^'^'^o  •=  "  ^1  ^^  ^ff>  *  <nJrJi>  ^  Y^  <kVT,-VT^>)  *   f  ,   (60) 

and  expresses  heat  production  by  the  fluctuations.   In  fact,  this 

Oh 
form    suggests  that  e  Is  the  turbulent  diffusion  flux  associated  with 

magnetic  field  gradients  (the  current),  and  <Siy^>  is  the  diffusion 

flux  associated  with  VT^,  namely,  the  turbulent  heat  conduction  term. 

This  interpretation  is  supported  by  the  result,   proven  in  Ref.   1^ 

[along  the  same  lines  which  lead  from  Eq.   (62)  to  Eq.   (65)  in  this 

paper],  that  <s^v^>  vanishes  when  <  -►  0  so  that  e  is  generated  by 
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motlon  due  to  resistive  instabilities,  and  by  the  result,  to  be  proven 
in  the  next  section,  that  e  vanishes  as  n  -*  0. 

Relation  (60)  can  be  used  to  guess  the  form  of  the  parallel 
component  of  e,  as  was  done  in  the  previous  section.  For  a  low 
pressure  plasma  J^  -  XB^  and  PqVT^  is  small.  Thus,  Eq.  (60)  implies 
/Ac-B^d^x  <  0.  Equation  (36)  remains  valid  here  as  well,  because 
redefinition  (55)  does  not  affect  the  parallel  component  of  e.  Using 
the  argument  given  in  Section  III  we  end  up  with  exactly  the  same  form 
as  in  the  incompressible  case,  namely,  Eq.   (38). 

Finally,  we  discuss  redefinitions  (55)  and  (56),  which  brought 

about   a   significant   simplification  of  the  evolution  equations 

(57)-(59).  The  effect  of  the  small  scales  enters  (in  addition  to  f, 

perhaps)   in  the  form  of  only  three  magnetic  surface  integrals  on  each 

surface;  two  associated  with  surface  components  of  e,  and  the  integral 
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or 


T  "^^l^ln^^^"  This  is  a  consequence  of  the  Grad-Hogan  theory' 
when  applied  to  our  equations.  The  small  diffusion  velocity  of  order  n 
is  now  Uq,  and  not  y^  which  may  be  cancelled  by  the  <PiV^>  term  in  Eq. 
(56).  This  term  acts  as  a  Gallilean  shift  in  the  frame  of  reference, 
and  the  redefined  e  is  related  to  <Vi'«B.|>  through  the  transformation 
governing  the  electric  field  when  observed  from  a  moving  frame. 

The  physical  meaning  of  u^  is  best  seen  within  the  context  of 
kinetic  theory.  Let  f(x,w,t)  be  the  distribution  function  of  the 
plasma,  where  w  is  the  microscopic  velocity  variable,  and  /fd-^w  -  n, 
the  number  density.  The  macroscopic  velocity  y  is  given  by 
nv  -  /wfd^w.  Let  us  write  f  as  the  sum  of  mean  and  fluctuating  parts, 


-26^ 


f  -  fjj  *  f ., .  Integrating  over  velocity  space  we  get  n  -  n^  ♦  n^  where, 
clearly,  n^  -  /f^d^w.  The  next  moment  yields 


IV  -  n^Up  +  J  wf^d^w  , 


"X  -  "n"o  *   wf.d^w  ,  (61) 


where,  by  definition,  u^  -  n^Vwr^d^w.  The  last  term  in  Eq.  (61)  has 
zero  mean,  and  we  have  n^UQ  -  <ny> .  This  coincides  with  the  definition 
of  Uq  in  Eq.  (56),  and  establishes  Uq  as  the  natural  macroscopic 
velocity  associated  with  the  mean  distribution  function  f„. 
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V.   TURBULENCE  WITH  ZERO  RESISTIVITY 

The  effect  of  the  turbulent  heat  conductivity  is  seen  most  clearly 
when  the  plasma  resistivity  vanishes.  This  is  the  case  we  now 
investigate.  The  same  scaling  as  in  the  previous  section  Is  assumed, 
where  now  the  small  parameter  is  attributable  either  to  k  (which  does 
not  appear  explicitly  in  the  following  calculations)  or  to  any  other 
cause  of  the  slow  diffusion  of  mean  fields.  We  will  show,  In  this 
case,  that  e  defined  by  Eq.   (55)  vanishes,  at  least  to  the  order  of 

l:^ilHlil- 

First,  we  consider  the  diffusion  of  a  scalar  function  i^  which  is 
advected  by  the  fluid, 


il  ■*■   W'V^^^   -  0  .  (62) 

3t 


Using  mass  conservation,  we  form  the  conservation  equation 
3(p4-')/9t  +  V«(p\J;v)  =  0.  Averaging  yields,  as  in  previous  such 
exercises, 


i-  (p^i^q)  *  V.(p^,^^u^  *   Po<'^iyi>)  -   0  ,  (63) 


3t 


where  <p^4;^v>  is  dropped  as  too  small.  In  Eq.  (63).  Uq  ^^  defined 
exactly  as  in  Eq.  (56).  We  now  notice  that  any  function  Fiii)  must 
also  satisfy  Eq.  (62),  as  well  as  an  equation  corresponding  to  (63). 
Writing  F(\|'o*<i'T)  •=  ^(^q)  *  F('i'o)'^-|  *  0{ii^),  where  the  dot  indicates  a 
derivative,  we  find  that,   up  to  an  error  of  order  \J)^  ,  F^  -  F(i^q), 
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F,  -  F(i|»  )4i^.  Thus,  Eq.   (63)  written  for  F  takes  the  form  (correct  to 
the  order  of  I'l'iY^  | ) , 

±   (PqF^)  *  V.[p^F^Uo  -  Po^(*o><*i:^1>^  -  0  •     ^^^^ 

Multiplying  Eq.  (63)  by  h^^)  and  comparing  it  with  Eq.  iSk) ,  we  get, 
after  using  Eq.   (58), 

<\1)^v^>'VFUq)   -  0  .  (65) 

Thus,  <'l'Ty^>  is  perpendicular  to  constant-iii^  surfaces.  We  now  identify 

il)  with  the  flux  function  of  B^  that  labels  its  toroidal  flux  surfaces, 
^o  -O 

Indeed,  the  average  of  V-(i^B)  -  0  yields,  to  order  1,  '•('J^qIo^  "  "^ 
which  implies  B^'Vii^  -  0.  In  all  the  following  discussion,  only  the 
component  of  <ij;^v^>  normal  to  \|)Q-surfaces  appears  (which  is  a 
consequence  of  the  Grad-Hogan  theory^°'^^ ) .  We  may  thus  view  Eq.  (65) 
as  implying  <'^iV^>  «=  0,  and  Eq.   (63)  reads 


3*0 

— -   ♦  u-'Vii)^  -  0  .  (66) 

3t    -°   o 


Consider  now  the  turbulent  diffusion  of  B,  which  obeys  the 
"frozen-in"  equation 


dB 

— 1  ♦  Vx(Bxv)  -  0  .  (67) 

at 
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The  field  b  ^  \pB  satisfies  the  same  equation,  and  its  average  (correct 
to  the  order  of  |y^Bj)  yields 


^  ('^o?o^  *  '''['J'o^?o''yo  -  £)]  -  0  .  (68) 


Using  Eqs.  (57)  and  (66),  we  find  £"^4)^  -  0.  Since  only  the 
components  of  t  tangential  to  t^^-surfaces  enter  the  equations,  we  may 
conclude  that  e  -  0  as  claimed. 

Returning  to  Eq.   (60),  we  find 


j  Po<s,y^>.VTQd3x  <  0  , 


(69) 


where  the  integral  is  taken  over  the  volume  of  the  plasma.  As  in  the 
previous  sections,  low-level  turbulence  with  negligible  correlation 
accross  flux  surfaces  implies 


?  ^^o 
<s,v.  >dS  =  -K"^  — -  (70) 

1  In  dV 


2 
on  every  mean  magnetic  surface,  where  K  is  a  non-negative  function. 

This  relation  suggests  that,  in  the  case  of  zero  resistivity,  turbulent 

heat  conductivity  acts  in  a  similar  way  to  classical  molecular  heat 

conductivity.    This  conclusion  may  have  to  be  modified  if  the 

resistivity  is  finite,  but  this  question  will  not  be  pursued  here. 
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VI .   CONCLUSION 

In  this  article  we  have  presented  a  dynamical  description  of  MHD 
plasmas  affected  by  low-level  turbulence.  Two-dimensional  fluid  flow 
was  also  discussed  as  a  simple  example.  While,  except  for  the  fluid 
case,  we  never  write  down  the  dynamical  equations  of  the  fluctuations, 
a  comparison  of  conservation  laws  for  the  full  variables  and  for  the 
mean  variables  makes  it  possible  to  derive  important  properties  of  the 
turbulent  diffusion  terms.  We  thus  find  that  turbulent  resistivity  is 

quite  different  from  its  molecular  counterpart,  and  can  explain 

q  ID 
spontaneous  toroidal  field  reversal  in  RFP  experiments.^'    Turbulent 

heat  conductivity,  on  the  other  hand,  is  more  similar  to  the  molecular 

effect  and  transports  heat  in  the  direction  opposite  to  the  mean 

temperature  gradient. 

An  issue  which  is  left  open  here  is  the  question  of  "cross 

effects";   the  effect  of  temperature  gradient   and   finite   heat 

conductivity   (which  unleashes  the  corresponding  instabilities)  on  the 

turbulent  resistivity,  and  the  effect  of  resistivity  on  the  turbulent 

heat  conductivity.  The  question  that  presents  itself  is  whether  there 

is  some  symmetry  relation  between  these  two  effects,  similar  to  the 

relation  which  exists  between  different  classical  diffusion  effects 

Oh 

according  to  Onsager's  theorem.  Our  expectation  is  that  the  answer  to 
this  question  is  in  the  affirmative,  and  base  our  opinion  on  the  fact 
the  Onsager  relations  apply  to  phenomenological  diffusion  laws, 
regardless  of  whether  their  origin  is  collisional  or  turbulent.  A 
credible  answer  to  this  question,  however,  has  to  be  based  on  the 
macroscopic  (fluid)  equations  of  the  turbulent  fluctuations. 
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FIGURE  1 


The  geometry  of  nested  streamlines 
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FIGURE  2 


Radial  profiles  of  (a),  the  toroidal  field  B^  and  (b),  the 
parallel  current  X  -  Jo*lo''l§ol^'  ^^^  shown  for  a^  -  1 , 
12.   Here  x  -  r/a   is  the  normalized  radial 


Bq  -  100,  E 


variable. 
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